Quantification of entanglement by means of convergent iterations 
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An iterative procedure is proposed for the calculation of the relative entropy of entanglement 
of a given bipartite quantum state. When this state turns out to be non-separable the algorithm 
provides the corresponding optimal entanglement witness measurement. 

PACS numbers: xx.xx 
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Entanglement is an important resource of quantum in- 
formation processing. Although there are quantum pro- 
tocols based on other features of quantum mechanics 
such as quantum superposition principle in celebrated 
Deutsch's searching algorithm GJ|, or protocols where the 
use of entanglement can be advantageous but is not es- 
sential, such as signaling through depolarizing channels 
with memory , most quantum protocols rely on the ex- 
istence of non-separable states. For practical purposes it 
is very important to quantify entanglement generated by 
realistic laboratory sources and thus evaluate the poten- 
tial usefulness of a given realistic source for the quantum 
processing/communication purposes. 

One of the measures of entanglement thoroughly stud- 
ied over the past decade is the relative entropy of entan- 
glement defined as j| 



E(a) = mi S(a\\p), 



where S is the quantum relative entropy, 



S(a\\p) = Tr(<7lner — a In/?), 



(1) 



(2) 



between states a and p; the infimum in Eq. (Q) is taken 
over the set of separable states. This functional is one 
possible generalization of the classical relative entropy 
between two probability distributions |I| to quantum the- 
ory. Let us mention that unlike in the case of entropy this 
generalization is by no means unique. This quantity can 
be given a geometrical interpretation as a quasi-distance 
between the state whose entanglement we are interested 
in and the convex set of separable states. E fulfills most 
of the requirements usually imposed on a good entan- 
glement measure and has other good properties. Most 
notably, on the set of pure states it coincides with the 
Von Neumann reduced entropy || ||, and is closely re- 
lated to some other measures of entanglement (7). Rel- 
ative entropy makes also a good entanglement measure 
for multipartite [g and infinite-dimensional || quantum 
systems. 

The analytical form of E is known only for some special 
sets of states of high symmetry jnj [HJ |l2|| . Generally 
one has to resort to numerical calculation. In a sense the 
problem resembles the reconstruction of quantum states 
using the maximum likelihood principle Jl3| , [l4| . Here the 
given input state a plays the role of experimental data; 
once this state is known, the statistics of any possible 
measurement performed on it is available. The solution 



can be obtained by means of several numerical methods. 
The formulation given in || is just an example corre- 
sponding to an implementation of the downhill simplex 
method. Its efficiency strongly depends on the dimen- 
sionality of the problem. 

In the procedure proposed here more analytical ap- 
proach will be adopted. We will derive a set of extremal 
equations for E and will show how to solve them by 
means of repeated convergent iterations. But this is not 
the only goal. The extremal equations indicate that there 
is a structure of quantum measurement associated with 
the extremal solution. The separable measurement ob- 
tained in this way specifies the extremal separable state 
and, significantly, it provides the optimal entanglement 
witness operator revealing the possible entanglement of 
the input state a. 

Let us denote p* the separable state having the small- 
est quantum relative entropy with respect to a. Let 
f(x,p*,p) — — x)p* + xp\ be the relative en- 

tropy of a state obtained by moving from p* towards 
some p. We are looking for the global maximum of a 
convex functional on the convex set of separable states. 
Two cases may arise. When a is separable the necessary 
and sufficient condition for the maximum of S is that its 
variations along the paths lying in the set of separable 
states vanish, 



dx ( ' 



p* , p) = 0, Vp separable. 



(3) 



When a is entangled S attains its true maximum out- 
side the set of separable states and we must carry on 
the maximization on the boundary. In that case Eq. (|^) 
holds only for variations along the boundary. It is well 
known that any separable state from the Hilbert space of 
dimension p — d ® d can be expressed as a convex sum 
of (at most) p 2 projectors on disentangled pure states 
(Caratheodory's theorem, see also §), 



p 

E 

k=l 



Vk)(<Pk\®\<Pk)(<fl\ 



(4) 



Here \<p\) and \<p\) are pure states (not normalized) of 
the systems 1 and 2, respectively. Now by taking squares 
of the projectors, \ipl' 2 )(ipl' 2 \ -> (\<pl' 2 ) 2 J) , one can 
remove the boundary and make condition (pf) universal. 
The derivation in Eq. (0) can easily be calculated using 



2 



an integral representation of the logarithm of a positive 
operator §. It reads 

g(0, p\ p) = Tr((p* + t)- l a(p* + t)- l Sp)dt 
= TrASp, 

where we denoted (1 — x)p* + xp = p* + Sp, and operator 
A has the following matrix elements in the eigenbasis 

{iA n )}of>* m, 



(X m \A\X n ) = 



log X n - log \ r , 

Xn A m 



■(A m |cr|A„). 



(6) 



Its meaning will be discussed later. The right hand side 
of Eq. (^) should vanish for all Sp preserving the form 
(||) of p* + Sp. There are p 2 such basic variations: tp\ ® 



ip'f, — ► {ip\ + 5ip\) ® (tpf. + Sip' 2 .). To make sure that the 
only constraint Tr(p* + Sp) = 1 is obeyed we will use a 
Lagrange multiplier A. 

Getting all things together we arrive at the extremal 
equations that read 



RI\<pI)(<pI\ = M<pI)(<pII 
RM)(<ti\ = \\<d)(<ti\, 



i. 



(7) 



where 



Rl = Tt a (A\<fi)(<? k \), 



Tr^K)^). (8) 



Bars denote projectors normalized to unity. Multiply- 
ing the first row of Eq. (Q) by |¥>|)(¥>||, the second by 
l < / , ifc)(Vfel) an d summing them separately over k we find 
that A = Tr^4p* = 1. Using this and modifying further 
the necessary condition we get the main formal result of 
this paper: The state having the smallest quantum rela- 
tive entropy with respect to a given state a satisfies the 
following 2p 2 equations, 



l\Rl = H){<pl\ 



k = i. 



(9) 



Unfortunately, solving such highly nonlinear operator 
equations by analytical means for anything but most triv- 
ial states seems to be out of question. One has to turn 
to numerics. We suggest to solve Eqs. (^) by repeated 
iterations starting from some randomly chosen separa- 
ble p. Let us note that the iterative procedure based on 
Eqs. @ belongs to the family of gradient-type algorithms 
of the form 4+ 1 = (dS{^) / dx^x^dS^) / dx k ) . Algo- 
rithms of this type are known to behave well; some of 
them were even proven to converge monotonically [fj"5| . 
They have found important applications in various opti- 
mizations and inverse problems. In our case we observed 
that the step generated by operators R k 2 m Eqs. @ was 
often too large — rather than converging to the station- 
ary point the algorithm would oscillate or diverge. If this 
happens the length of the step can be made smaller by 

12. 

mixing the operators R k with the unity operator: 



R 



1,2 



(1 



|«^' 2 )/(1 



(10) 



Indeed, when a is sufficiently small the algorithm con- 
verges monotonically. This can be seen by considering 
an infinitesimal step with a«l. It is convenient to split 
one iteration of Eqs. (||) into two subsequent steps corre- 
sponding to the two rows of Eqs. (|J) (projectors of only 
one of the subsystems are updated at a time). The two 
steps are completely symmetrical, so we will consider an 
infinitesimal iteration on, say, the projectors IVfc'XVfc'l 
of the first system obtained after the i-th iteration. Wc 
want to show that after one such step the quantum rel- 
ative entro py is never increased, S(a\\p 1+1 ) < S , (cr||p I ). 
Using Eq. ( |l0| ) in Eqs. (^|) we get to the first order in a, 



= (1 - a)^ 1 + ap, 



(11) 



where p = \ Z k (*%\<P%) (^l+K 1 ) (viW)®lv* > (vtl 
and thus 

S(a\\p i+1 ) - SHI,/) cx = ! _ ^i- (12) 

oa 

It remains to show that TrA'p > 1. Let us denote 
K = (fk\fk)(<Pk\<pf)- Notice that J2 k X[ = 1 by the 
normalization of p 1 . Then by using the Swartz inequality 
and the concavity of the square function we obtain 



Tr J 4 i p = ^Aj fe Tr(^ i |^}(^|^ i ) 



> 



L[i*(j#R>WD 



k 



(13) 



= (TrAV 1 ) = 1, 

which completes our proof. This means that with a suf- 
ficient amount of regularization (|TFJ) the algorithm (j^) 
converges monotonically. In practice, the parameter a 
need not be very small. In 2 (8> 2 and 4 ® 4 dimensional 
problems we tried monotonic convergence was observed 
even with a of the order of unity. 

Now let us go back to the extremal equations @ . The 
left hand sides are generated by the operator A, which de- 
pends on a through Eq. Let us assume for now that a 
is an entangled state. Then A represents the gradient of 
the quantum relative entropy <S'(o'|| j o) at p* — the separa- 
ble state closest to a. Loosely speaking, the states giving 
the same expectation, TrA(a)p = const., form hyper- 
planes that are perpendicular to the line connecting a 
and p*. Since TrA(a)p* = 1 and p* lies at the boundary 
of the set of separable states, the conjecture is that the 
operator A is up to a shift of its spectrum a witness op- 
erator |l6|, [ItJ detecting the entanglement of a. In the 
following we will show that this is indeed the case, and 
that the operator 



W(a) = 1 - A(a) 



(14) 



is indeed the optimal witness of the entanglement of a. 
The mutual relationship of a, p* , and the states detected 
by W is shown in Fig. |l|. 
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FIG. 1: Mutual relationship between an entangled state a, 
the separable state p* closest to it in the sense of quantum 
relative entropy, and the entangled states detected by the wit- 
ness operator W. 



First we will show that TiAp < 1 if p is separable. To 
this end let us note that 



1 - TrAp 



df 
d.r 



lim 

x—>0 



(0,P*,P) 

S(a\\(l-x)p*+xp)-S(a\\p*) 



(15) 



Now, since both p and p* are separable states, so is 
their convex combination (1 — x)p* + xp. But p* mini- 
mizes S'(o'llp) over the set of separable states. Therefore, 
S(a\\(l-x)p* +xp) -S(a\\p*) > 1. This holds for all x 
so we have, 

TrW(cr)p = 1 — TrA(a)p > 0, Vpseparable (16) 

This already means that W is an entanglement witness 
operator. To show that W detects a we will again make 
use of Eq. (|l5|) with p now being substituted by the en- 
tangled state a. Now, because of convexity of S, 



S(a\\(l-x)p* 



S(a\\p*) 



< S(a\\p*) < 0. 

(17) 

The last inequality follows from the assumed non- 
separability of a. Eq. ( |l7| ) also holds for any x so we 
obtain 

TiW(a)cr = 1 - TrA(o-)cr < 0, Vcr entangled (18) 

which we set out to prove. 

Possible applications of our algorithm are twofold: 
First, it can be used for checking whether a given state 
is separable or not. Second, it can be used for quantify- 
ing the amount of entanglement the state contains. As 
a test of separability we tested the algorithm on many 
randomly generated separable and NPT states of dimen- 
sions 2 (g> 2 and 4 ® 4; typical results are summarized in 

Fig. |. 
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FIG. 2: Test of separability. It is shown how the calculated 
relative entropy of entanglement of several randomly gener- 
ated separable states approaches zero in the course of iterat- 
ing. The ordinate is labeled by the precision in decimal digits. 
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FIG. 3: The accuracy of the calculated relative entropy of 
entanglement is shown after a given number of iterations for 
six Werner states with / G [—0.05, —0.8]. The ordinate is 
labeled by the precision in decimal digits. 



Recently, another numerical test of separability has 
been proposed 1 18 consisting of a hierarchy of gradually 
more and more complex separability criteria that can be 
formulated as separate problems of the linear optimiza- 
tion theory. The algorithm we propose is much more 
simple. There is just one set of equations to be solved 
by repeated iterations and after that one finds not only 
whether the input state is entangled but also how much. 



Unfortunately, explicit formulas for the quantum rel- 
ative entropy are known only in very few cases. One 
such exception is the family of Werner states defined as 



4 



follows, 



Pw = 



d-f fd-l 

-1 + -rr^. ttF, 



d(d 2 - 1) d{d 2 - 1) 



(19) 



where F is the flip operator <8> ip 2 ) = ip 1 ® ip 2 and 

/ = Trp w F G [—1, 1] is a parameter. Fig. || shows the per- 
formance of our algorithm for several entangled Werner 
states of dimension 4 (g) 4. It is worth mentioning that 
the optimal entanglement witness W for the detection 
of Werner states generated by the operator A Eq. (||) 
is simply W = 1 — 2|^_)($_|, where is the singlet 
state. The expectation value of W is then a renormalizcd 
singlet fraction. 

The curves appearing in Figs. || and || suggest that 
the convergence of the proposed algorithm is faster than 
polynomial but slower than exponential. In some cases 
such as that shown in Fig. || the convergence is nearly 
exponentially fast. This is just a qualitative statement 



since we did not attempt to do any optimization of the 
length of the iteration step. In all probability such op- 
timization would result in further speedup compared to 
our examples given in Figs. ^| and |^. 

In conclusion, we derived a convergent iterative algo- 
rithm for the calculation of the relative entropy of en- 
tanglement. It can be used for checking whether a given 
input state is entangled. If it is, the algorithm calcu- 
lates its relative entropy of entanglement, finds its closest 
separable state, and provides the optimal entanglement 
witness measurement. 
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